Formation of H3 by radiative association of H2 and H in the interstellar medium 
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We develop the theory of radiative association of an atom and a diatomic molecule within a close- 
coupling framework. We apply it to the formation of HJ after the low energy collision (below 0.5 eV) 
of H2 with H~. Using recently obtained potential energy and permanent dipole moment surfaces 
of H;^, we calculate the lowest rovibrational levels of the H;^ electronic ground state, and the cross 
section for the formation of by radiative association between H[~ and ortho- and para-H2. We 
discuss the possibility for the Yi^ ion to be formed and observed in the cold and dense interstellar 
medium in an environment with a high ionization rate. Such an observation would be a probe for 
the presence of H~ in the interstellar medium. 



I. INTRODUCTION 

Many chemical reactions in the interstellar medium 
(ISM) are powered by cosmic rays: atoms and molecules 
(mainly molecular hydrogen) are ionized by radiation 
that provides sufficient energy to initiate a chain of chem- 
ical reactions in interstellar clouds leading to the synthe- 
sis of polyatomic molecules. About fourteen positive ions 
have been observed and identified in the ISM, including 
the simplest triatomic one, H^. This ion plays an impor- 
tant role in chemistry and evolution of interstellar clouds 
[U, 01, as its abundance is strongly related to the pro- 
duction of in the ISM. In contrast, the existence of 
stable negative ions has long been thought impossible due 
to the presence of UV radiation in the ISM. In particu- 
lar, the ion, which has a single bound state, has not 
been directly detected yet so that its presence, though 
probable, is still controversial Even if it has been 
proposed some time ago the surprise came with the 
detection of negative molecular ions in dense molecular 
clouds or carbon star envelopes. There are several car- 
bon chain ions that have been detected so far: CeH^ 
[113, C4H- i, CgH- [p, C3N- [g, C5N- gg, and 
CN~ [ll| . All these species are closed-shell systems with 
a quite large electron binding energy (3.8 eV for CgH"), 
and a large permanent dipole moment of several debye. 
Although it is generally accepted that in the ISM these 
ions are formed after radiative capture of an electron by 
the parent neutral molecules (12] . no detailed quantum 
mechanical treatment has been developed yet for model- 
ing the process. 

As the negatively charged counterpart of H;^ , the 
ion is predicted to be stable by about 0.013 eV [isl. IT^. 
It has never been detected so far in the cold regions of 
the ISM, while Wang et al. observed it by mass 

spectrometry of laboratory plasmas. The H;^ ion cannot 
be formed by radiative attachment from the unstable II3 
molecule, so that one has to consider alternative paths: 
the three-body recombination (TBR) H2 + + X ^■ 
-I- X , or the radiative association (RA) II2 -I- — > 
-|- hu) . The TBR mechanism must be the dominant 
process of Kg formation in the laboratory plasma but it 
is most probably inefficient in the ISM because of low 



H2 densities. The goal of this study is to investigate the 
possibility to form H;^ by the radiative association of II2 
and H" in low temperature (<150 K) environments. 

The collisions between H~ and H2 have been studied 
both theoretically [l^, [l^ and experimentally [itI . IT8| . 
only for collision energy larger than 0.5 eV. In contrast, 
the structure of the Hg ion has been rarely explored in 
the past [1^ . Recently, we have calculated (see Ref. |14| , 
hereafter referred to as paper I) a new accurate poten- 
tial energy surface (PES) for the Hg" electronic ground 
state, whose accuracy was improved compared to the 
previous ab initio calculations ^13 , I61] . In addition, we 
have determined for the first time the Hg permanent 
dipole moment surface (PDMS), which is needed for the 
RA calculations. The ion is well represented as a 
loosely-bound H2 • • • H~ complex with several rovibra- 
tional states, bound at most by about 70 cm~^. We have 
also found that there are a number of predissociation res- 
onances, which can be described as excited rovibrational 
states (j, v) of II2 perturbed by H" , coupled to the dis- 
sociation continuum Il2(j', f')+H^ with energy of the 
dimer state E{j',v') lower than E{j,v). 

The paper is organized as follows. In the next section 
we discuss the geometry of the Hg" molecule. In section 
mil we introduce the rovibrational wave functions of the 
molecule for bound and continuum states. The theory of 
radiative association of a dimer and an atom is developed 
in section IIVI Finally, in section |V] we present results 
of numerical calculations of the RA cross-section and of 
the rate coefficient for H;^ formation, and we discuss the 
possibility to observe the H^ ion in the ISM. 



II. REPRESENTATION OF THE H2 
COMPLEX AT LOW ENERGIES 



H" 



The Hg" molecule is composed of three identical nuclei, 
described in principle within the CNPI (complete nuclear 
permutation inversion) group D^h [T^ . However, due to 
the loosely-bound nature of the ion in its electronic 
ground state |13,, 14] , the exchange probability of the H^ 
proton with the dimer protons is negligible. Therefore, 
the low-energy collision between H2 and H~ can be stud- 
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ied as an inelastic collision of two structured particles. 
We first define two coordinate systems: The space-fixed 
(SF) frame with axes {x,y,z), and the body-fixed (BF) 
frame with axes {X, Y, Z) whose orientation in the SF 
frame is given by the Euler angles ae, fie and 7e. The 
Z-axis connects the center of H2 molecular axis with the 
nucleus of H~, and its X-axis is in the plane of the three 
nuclei. The orientation of the dimer with respect to Z is 
given by the azimuthal angle 6. The natural coordinate 
system (CS) associated to the BF frame is the Jacobi 
CS (i?, r, 9) where R is the distance between H2 and H~ 
along the Z axis, and r the internuclear distance of H2. 

In the framework of this " super-dimer" approximation, 
the quantum numbers of the H2 dimer are treated as good 
quantum numbers. The quantum states of the H2 • • • H~ 
complex at low energies are characterized by the quan- 
tum numbers (J, M, j, u, il^): J and j, associated to the 
total angular momentum J of the complex and to the 
angular momentum j of the dimer, the vibrational state 
V of H2, the absolute values M and Q. of the projection 
of J on the SF z axis, and of j on the BF Z axis, and 
the intrinsic parity ± of the wave function with respect to 
the reflection cr^ through the plane containing the Z-axis. 
For 7^ the $7+ and states are degenerate in this 
approximation. If the II2 ■ • • complex is bound, an ad- 
ditional quantum number vt characterizes the vibration 
along the Z axis. One defines the basis set | J, u, j, Vl) 

\j,v,3,n)= (1) 



2J+1 
An 



or, in the equivalent form 



\J,v,j,n) 



(2) 



2J+1 



87r2 



Dij^iae, Pe, 7e)] * e5^(cos 0)xvj{r) 



e5^(cos( 



where the normalized function 

Pj^ {cos 9) / ^/Att [i^l is proportional to the associ- 
ated Legendre polynomial and describes the rotational 
state of the dimer, and Xvji''') is the vibrational wave- 
function of II2. The Wigner function Df^Q{ae, Pe,7e) 
is associated to the rotation of the complex in the SF 
frame. The function above is not yet symmetrized with 
respect to and, correspondingly, ft can be positive or 
negative. Assuming that the PES V{R, r, 9) is known 
(see paper I), one derives the interaction matrix V 
(diagonal with respect to 51) written in the above basis 
(after integration over r and 0, denoted by the subscripts 
of the angled brackets) 



J 

vjfl,v'j'fl 



,{R) = {Jvjn\V{R,r,9)\Jv'fn)r^0dnn' • (3) 



As i? — > cx), the diagonal (in v and j) elements of this ma- 
trix become equal to the energies e^j of the rovibrational 
quantum states \vj) of H2. The diagonalization of the V 
matrix at every R value yields a set of adiabatic poten- 
tial energy curves, correlated at large R to the sum of e^j 
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FIG. 1. (Color online) The adiabatic curves of H;^ as func- 
tions of the Jacobi coordinate R, and the scheme for radiative 
association at a low coUisional energy E. Horizontal lines in- 
dicate positions of several bound states (reported in Table H]) 
and resonances. 



and of the ground state (Fig. [T]). These curves can 
accommodate bound levels corresponding to the quanti- 
zation of the vibrational motion along the R coordinate. 
The H2 nuclear spin = (para-H2) 01 id — ^ (ortho- 
H2) determines the possible values of j, namely even in 
the former case, and odd in the latter case. The energy 
of the state does not depend on the quantum number M 
which will be omitted in the following. 

The wave functions of Eq. 1^ must be properly sym- 
metrized with respect to and total inversion E* , cor- 
responding to the symmetry operators of the Coov group 
relevant to the present approximation. The |J, f2) 
functions obey the transformations 



ay\j,v,j,n) = \j,v,j,-n) , 
E*\j,v,j,n)^i-iy\J,v,j,~n), 



(4) 



so that we define the appropriate symmetrized combina- 
tions 

\j,v,j,n^) = i\j,v,j,n,)±\j,v,j,-n))/V2. (5) 

From these equations, the total parity of the -1- and — 
states is given by the value of (—1)'^ and (—1)''+^, re- 
spectively (also referred to as the e and / parity for linear 
molecules 21]). The permutation (12) of the the two nu- 
clei of II2 requires another symmetrization of the dimer 
wave functions 



{i2)\j,v,j,n) = {-iy\j,v,j,n), 



(6) 



so that we can write the following four wave functions 
labeled by the irreducible representations (irreps) of the 
C2v group r = ^1,^2,^1,52 (see Eq. (94.18) of 0) 
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(with n>0): 

\A,,j,vj,n}=M{\j,j,n} + {~iy\j,v,j,n} 
+i-iy\j,j,-n) + i-iy+^\j,v,j,^n)) (7) 

\A2,j,v, J, n) = AA(| J, J, 0) + (-1)^ I J, v,j, n) 
-(-1)^1 J, J, -17) - i-iy+^\j,v,j,-n)) (8) 

\B,,j,v,j,n) =j^{\j,j,n) - {-iy\j,v,j,n) 
-{~iy\j,j,-n) + {-iy+^\j,v,j,-n)) (9) 

\B2,J,v,j,n) - {^iy\J,v,j,n) 

+{-iy\j,j-n) ~ (-i)-^+^|j,«,j, -r!)) . (lo) 

The normahzation factor M above is 1/4 or l/(2-\/2) for 
= and 17 > respectively. Some functions above 
are identically equal to for certain combinations of J, j, 
and 17, which means that the corresponding irreps are 
not allowed to these combinations. 



III. BOUND AND CONTINUUM WAVE 
FUNCTIONS OF THE Ha H COMPLEX AT 
LOW ENERGIES 

The super-dimer approximation above is useful for 
defining basis functions with the appropriate symmetry. 
However, the quantum numbers of the H2 internal state 
can be used at large separation, but not when H2 ap- 
proaches H^. An accurate description of the wave func- 
tions |r, J, E) of the H2 • • • H~ complex with energy E 
(E being either the collision energy, or the bound state 
energy, of the complex) is obtained within the close- 
coupling framework for which we recall below the main 
steps. The close-coupled expansion on iVtot symmetrized 
channel functions |r, J, a) is written as 



iVtot 



V,J-E) ^\^Ca{E)\T,J,a-E). 



(11) 



where a is the channel index. The channel func- 
tion |r, J, a; E) are the independent solutions of the 
Schrodinger equation with appropriate boundary condi- 
tions. They are written as a linear combination of iVtot 
basis functions |r, J, u, j, 17^) 

Wtot 

|r, J, a;E) = Y. ^kAR; E)\T, J, Vk,jkMk) ■ (12) 

k=l 

In these expressions, the index a can be assigned to the 
asymptotic limit of the channel wave function in the BF 
frame at large distances, i.e. {a = Va,ja,^a)- The num- 
ber of channels iVtot is determined by convergence check 
on the calculated physical property, namely the binding 
energies or the RA cross sections. 

The Fk,a expansion coefficients are solution of the set of 
coupled differential equations written (in a matrix form) 
in the BF frame as 



I+{e-El)+ W{R) 



F{R)^0, (13) 



where m is the reduced mass of the complex, and e the 
diagonal matrix of the H2 energies eyj . The interaction 
matrix W is the sum of the interaction potential ma- 
trix V (Eq. ([3]) and of the matrix R (diagonal in v and 
j) describing the atom-diatom relative rotation with the 
orbital angular momentum operator Z = J — j. 



r,j,v,j,n\P\T,j,v',j',n') 



^vj.v' j' 



(14) 



A detailed expression of these matrix elements in the BF 
frame can be found in Ref. 22 1. The operator 1? couples 
channels with different values of 17, through long-range 
terms which dominate the potential energy terms due to 
their character. Therefore in order to treat conve- 
niently a scattering process, it is usual to choose the rep- 
resentation in the SF frame where the scattering channels 
are fully decoupled, and where the P operator is diagonal 
(with eigenvalues £{£ + 1)): 



^vj£,v' j'£' 



1) 



2mi?2 



'^vj£.v'j'£' 



(15) 



Similarly to Ea. ([TT|) . the total wave function \T, J;E) in 
the SF frame is expressed as 



\r,J;E) = J2Co.iE)\r,J,a;E), (16) 



a=l 



where the channel wave functions jF, J, a; E') are ex- 
panded on the basis set jF, J, v,j, i) 

Ntat 

\T,J,a;E) ^Y.P-AR-^^WA,v.,j.J.) ■ (17) 

Like the channel index a in the BF frame, the index a can 
be assigned at large distances to the asymptotic channel 
i.e. {a = Va,ja,£a)- The F^.q coefficients are solutions 
of Eq. ([T3|) where the matrices are expressed in the SF 
basis set. The matrix elements of V in the SF frame are 
related to those in the BF frame according to 

J 

n=o 

where (i'|17)'"'' is a matrix element of the transformation 
between BF and SF angular basis sets 



(£|17) 



rj 



i-iy+''V2£ + TA2 - 5n,o 



j J £ 
17 -17 



(19) 

The total parity (—1)'^ and (—1)'^+^ for -|- (or e) and — 
(or /) states imposes that the sum in Ea. ([T71) is restricted 
to the j and £ values satisfying (—1)'' = (— 1)-'+^ and 
(-l)"'+i = (-l)^'+^ respectively. 

At this step, the coefficients of channel superposition 
in Eq. (in]) are still unspecified. They are determined 
by applying appropriate boundary conditions when the 
set of differential equations of Eq. (IT51) is solved. Bound 
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states are characterized by imposing the standard bound- 
ary conditions for each channel a 

F,AO;Et) -0 V«e[i,iVtot] 

F^MEt) VkG [l,7Vtot] . ^ ^ 

These quantization conditions are fulfined only for dis- 
crete energies Et , corresponding to a discrete vibrational 
level Vt of H;^ characterizing the number of nodes (in the 
R coordinate) of the radial wave function in the dominant 
channel o'q for given F and J. We used the renormalized 
Numerov method [1^ to solve the coupled equations of 
Eq.([T3]) written in the SF frame, so that the resulting 
superposition coefficients C„* ensure that the outward- 
propagated and inward-propagated wave functions are 
equal with identical derivatives at a well-chosen match- 
ing distance [23*]. Equations (|16p and (ITTI) are recast for 
the wave function of an bound level |F, J, Vt) as 

Ntat 

\r,j;a'„vt) = $r&(^)|r, (21) 

k' = 1 

Ntat 

= E C:rF,,,^,{R;Et). (22) 

q' = 1 

As an illustration, the calculated energies of lowest 
bound levels for J = to 6 corresponding to the binding 
of para-H2(f = 0, j = 0) with H~ are displayed in Ta- 
ble [H labeled with vt and with the approximate quantum 
numbers v,j,£ of the corresponding dissociation thresh- 
old. This is justified by the weak interaction of the lowest 
channel {v — 0,j = 0,i) with the other ones, as already 
noticed in paper L There is a single dominant compo- 
nent in the expansion of Eg. (1^1]) . The present energies 
for J = are found in agreement within 0.5 cm~^ with 
the ones obtained in paper I using a different integra- 
tion method. Beyond J = 6 no bound level can exist 
anymore below the dissociation threshold. Bound lev- 
els with energies below the ortho-H2(f = 0,j ~ 1)+H~ 
threshold are reported in Table |lll In the C2v frame- 
work above, these levels correspond to the binding of 
a para-(ortho)-H2 with H~, and all other bound levels, 
which are found at energies below thresholds with non- 
zero even (odd) values of j correspond to predissocia- 
tion resonances, which will decay into the continua of 
the }i2{v = 0,j')-|-H~ thresholds with j'even (odd) and 
0(1) < j' < j. All these levels are included in the RA 
cross section calculations below. 

Note that such values could be used for the search for 
H;^ rotational transition lines in the ISM absorption spec- 
tra in the mm-wavelength range. As the component of 
the permanent electric dipole moment of along the R 
axis exceeds by two orders of magnitude the transverse 
components (paper I), we can follow the notations of 
Ref.[21j for electric dipole transitions in a linear molecule. 
Q lines (A J = 0) will occur between e and / states, and 
will not be observed between levels of Table H] P and R 
lines (AJ = ±1) will connect levels with the same parity 
e or /. 



J,v,j,£,vt 


Energy 


J,v,j,£,vt 


Energy 


0,0,0,0,0 


-71.2 


1,0,0,1,1 


-23.8 


1,0,0,1,0 


-68.2 


2,0,0,2,1 


-20.2 


2,0,0,2,0 


-62.4 


3,0,0,3, 1 


-15.0 


3,0,0,3,0 


-53.8 


4,0,0,4,1 


-8.6 


4,0,0,4,0 


-42.8 


5,0,0,5,1 


-1.7 


5,0,0,5,0 


-29.7 


0,0,0,0,2 


-5.4 


6,0,0,6,0 


-15.3 


1,0,0,1,2 


-4.5 


0,0,0,0, 1 


-25.7 


2,0,0,2,2 


-2.7 






3,0,0,3,2 


-0.4 



TABLE I. Computed binding energies (in cm~^) of the vibra- 
tional levels Vt ordered as rotational progressions in J, of the 
para-H2-H~ complex with respect to the lowest dissociation 
limit H2(i' = 0, j = 0)+H~. Levels are labeled with the ap- 
proximate quantum numbers f = 0, j = 0, and £ = J. In the 
super-dimer picture, they are all of e parity. 



J,v,j, 




Vt 


Energy 




i, 


Vt 


Energy 


0,0,1 


1 





-152.7 


0,0,1 


1 


2 


-33.2 


1,0,1 








-149.3 


1,0,1 





2 


-31.6 


2,0,1 


1 





-142.6 


2,0,1 


1 


2 


-28.4 


3,0,1 


2 





-132.6 


3,0,1 


2 


2 


-23.9 


4,0,1 


3 





-119.7 


4,0,1 


3 


2 


-18.2 


5,0,1 


4 





-104.0 


5,0,1 


4 


2 


-11.8 


6,0,1 


5 





-85.8 


6,0,1 


5 


2 


-5.0 


7,0,1 


6 





-65.6 


1,0,1 


2 





-13.3 


8,0,1 


7 





-44.0 


2,0,1 


3 





-10.0 


9,0,1 


8 





-21.7 


3,0,1 


4 





-5.3 


0,0,1 


1 


1 


-77.8 


0,0,1 


1 


3 


-10.7 


1,0,1 





1 


-75.3 


1,0,1 





3 


-9.8 


2,0,1 


1 


1 


-70.5 


2,0,1 


1 


3 


-8.1 


3,0,1 


2 


1 


-63.3 


3,0,1 


2 


3 


-5.6 


4,0,1 


3 


1 


-54.2 


4,0,1 


3 


3 


-2.7 


5,0,1 


4 


1 


-43.2 


0,0,1 


1 


4 


-1.7 


6,0,1 


5 


1 


-31.0 


1,0,1 





4 


-1.3 


7,0,1 


6 


1 


-17.9 


2,0,1 


1 


4 


-0.7 


8,0,1 


7 


1 


-4.9 


1,0,1 





5 


-0.07 


1,0,1 


1 





-13.4 


4,0,1 


4 





-0.5 


2,0,1 


2 





-10.2 


1,0,1 


1 


1 


-0.1 


3,0,1 


3 





-5.8 











TABLE II. Computed binding energies (in cm~^) of the vi- 
brational levels Vt of the ortho-H2-H~ complex ordered as 
rotational progressions in J, with respect to the lowest disso- 
ciation limit H2(« = 0, j — 1)4-11"). Levels are labeled with 
the approximate quantum numbers ii = 0, ji = 1, and £. In the 
super-dimer picture, levels in the upper part are of e parity, 
and in the lower part of / parity. 



In a scattering process with energy E, the number of 
energetically open channels No gives the number of phys- 
ical solutions of Ea. (fT3|) . The remaining Nc = Ntot — A^o 
closed channels should lead to non-physical long-range 
behavior since the related wave functions diverge asymp- 
totically. Accordingly, the sum in Eq. is restricted 
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to the No open channels 

\T,J,a;E)^ (23) 



a=l 



K=l 



and the usual boundary conditions are applied for each 
a to the set of coupled equations 

i^^«,„(0;£;) =0 VaG[l,iVtot] 



E) 



R^oc 



R—>oo 



> Vae [No + l,Ntot] 

> j{R, k^)5^^a - n{R, ki^)K^,^a{E) 

Vae [l,No] . 

(24) 

The standard collision matrix K(i?) summarizes the 
channel interactions. The functions j{R,k^) and 
n{R,kK,) are the regular and irregular spherical Bessel 
functions with momentum fc^ with respect to the 
H2(wk, jK,^K)+H~ limit. To be fully determined, the to- 
tal wave function \T,J,E) must also specify the initial 
state of the scattering process, i.e. the entrance channel 
ao, through standard outgoing wave normalization 



C^o{E) = [I-^K{E)]-^^^ 



(25) 



IV. 



CROSS-SECTION FOR DIMER-ATOM 
RADIATIVE ASSOCIATION 



The RA of H~ with H2 in an initial rovibrational level 
(v,j), creating Hj^ in a vibrational level vt, with respect 
to a dissociation limit H^+H2(w', /) 



H 



u^ir',j',vtiv',j')) + hut, (26) 



is schematically shown in Fig. ^ioi v — 0,j — 0. A 
similar picture holds for H~ and ortho-H2(w = 0,j = 
1) association. We limit our study to these two states, 
which are the only ones populated in the cold and dense 
molecular ISM (see, for instance, Ref. jl^])- The two 
partners approach each other with a collision energy E 
above the threshold. A photon with energy hujt — E+Ey^ 
is emitted to stabilize in the vibrational level vt with 
binding energy Et, with respect to the H2(w — 0,j = 
0)+H- or the H2(u = 0,j = 1)+H- limits. 

In order to obtain the RA cross-section and rate coef- 
ficient, we modify the theoretical approach developed by 
Herzberg [2^, which was later used by several authors 
[26'-'28'| for RA in diatomic molecules, and for photo- 
association of cold atoms [1^ H^]. This formalism is 
somewhat similar to the treatment one of us used pre- 
viously for the photodissociation of van der Waals sys - 
tems, which can be considered as the inverse process [31| . 
The initial and final states are specified in the SF frame 
by \T, J, E,ao) and |r, J, iij, ag). The entrance channel 
index ao is correlated to the asymptotic quantum num- 
bers wq, jo of the initial H2 level, and £0 specifying the 
initial coUisional state. The final state of the created 



bound molecule has a multichannel nature, but it 
is labeled for convenience with the vibrational index Vt, 
and the index u'q of the dominant channel in its expan- 
sion (Eq. ipT]) '). also related to the asymptotic quantum 
numbers a'g = fo,jo'^o (^^ Tables HI and HH) . 

The Einstein coefficient A^/,^^"^ {E, Vf) for the process 
of Eq. ((26|) is written (in atomic units) p8l | 



3c3 



E 

CT=0,±1 



{r'J'vta'o\ti''\TJEao] 



27) 



where fj,'^ = ijP,ij,^^ are the three components of the per- 
manent dipole moment operator fi in the SF frame. The 
/i"^ components in the BE frame are related to the ones 
in the SF frame (determined in paper I) by the Wigner 
rotation functions D depending on the Euler angles ae 
and l3f. 



J2 [Ci,A(ae,/3e,0)]V' 



(28) 



A=0,±1 



We obtain 



a'r,;ao ^ ' ^' 



Au4 
3c3 



(2J+ 1)(2J' + 1) X 



E 

IT=0,±1 



J 1 J' 



{T\J'■,vt,a'a\^l'^\T,J■E, 



(29) 



A=0,±1 

which becomes, after averaging over the M values 



(2J' + 1) X 



Y {T',J'■,Vt,a'^\^i^\T,J■E,aQ) 



A=0,±1 



(30) 



We now introduce the expressions for the total wave func- 
tions (Eqs. dm [231)) 



No 

E CZ^'{E)dXf-^\E,vua', 



A=0,±1 
2 



a=l 



with 



^Aa {E.vt.a^) = 
E E(*rfc(^)lM^.«(i?)I^^K.a(i?;£;))«. 



(31) 



(32) 



1 = 1 K=l 



The subscript R at the angled brackets denotes the inte- 
gration over the R coordinate (33 | . The matrix elements 



of the i?-dependent transition dipole moment in the SF 
basis are 

fi^,,jR) = (33) 

O'=0 0=0 

related to the matrix elements in the BF basis 



J 1 J' 



{Xv'Ar)Qr'mfi\R,r,e)\x.,{r)ef{0))r,e. (34) 

The subscripts r, 6 at the angled brackets denote the in- 
tegration over the r and 9 variables. Note that the ini- 
tial and final wave functions are not necessarily described 
with the same number of channels, i.e. A'tot and Nl^^ can 
be different because J 7^ J and F 7^ F'. We use energy 
normalization for the initial continuum state. 

The probability P^,'^^'^''' {E,vt) of an RA event start- 



ing in the |F, J; E, ao) state with momentum k = \j2mE 
toward the final state |F', J'; vt, ap) is given by the Ein- 
stein coefficient divided with the current density of inci- 
dent particles, which is l/(27r) for the energy normalized 
wave function. The corresponding RA cross-section is 
then expressed as 



(35) 



In the following, the indexes ag and ctg are assigned 
to their asymptotic labeling (fo,jo,^o) and (fo,Jo,^o). 
Physically, the entrance channel is determined asymptot- 
ically by the initial rovibrational H2 level (wq, Jo); ^^at 
a sum over must be performed. Since the nuclear spin 
is conserved during the RA process, there is no need to 
include the nuclear spin degeneracy factor. All possible 
final states should also be included for the computation 
of the total RA cross section. For each initial J value, the 
contributions from each final value J' = J, J ± 1 must be 
added together. At a given collision energy E^ the sum- 
mation over all possible values of J must be performed 
as well. Finally, we are interested in the formation of 
in any of its stable bound levels Vt-, so that the total cross 
section is obtained as follows 



E 



r'j',rj 



As reported in paper I, the component /i°(-R, r, 9) (along 
Z-axis) of the H;3" dipole moment in Eq. is by two 
orders of magnitude larger than the two others compo- 
nents ii^^{R,r,9). Therefore, the contribution from the 
transverse components is neglected. The total cross sec- 
tions croo{E) and aoi{E) for the RA of H~ with para-H2 
{v = 0,j = 0) and with ortho-H2 {v = 0,j = 1) re- 
spectively, are shown in Fig. We already mentioned 
that at low collision energies (below 400 cm~^), a sin- 
gle component dominates the multichannel expansion of 




o 10 



Energy (cm ) 



FIG. 2. The RA cross-section (in atomic units) starting from 
para-H2 {vq = 0, jo = 0) (full line), and from ortho-H2 {vq = 
0,jo ~ 1) (dashed line) as a function of collision energy E 
above the corresponding thresholds. The inset enlarges the 
region of Feshbach resonances induced by the bound states of 
the closed channel para-H2(i' = 0,j = 2)-|-H~ (see Fig. [T]). 



the initial and final wave functions. Therefore, the ap- 
proximate selection rules Vi ^ Vf = Vi,ji — ?► j/ = ji, 
and J — J' = J ± 1 hold for the RA process (with the 
appropriate selection rule for the parity). As expected 
from Eq. (PT|) . the cross section for RA with ortho-H2 is 
larger than the one with para-H2, as the former species 
has a deeper potential well, and has more bound states 
than the latter. The bumps near 1.5 cm^^ visible in 
both curves are due to the enhancement of the prob- 
ability density at the top of the centrifugal barrier for 
low J. Moreover, one shape resonance is predicted near 
4.3 cm~^ in aoo{E) and near 6.5 cm~^ in aQi{E) associ- 
ated with non-zero J values. Once the collision energy 
reaches 285 cm"-'^ at which the j = 2 rotational state 
of para-H2 can be populated, a series of Feshbach reso- 
nances induced by the bound states of the closed channel 
para-H2(w = 0, j = 2)+H~ appears in the cross section. 

The RA rate coefficient /cra(2^) is obtained by a stan- 
dard integration over the Maxwell-Boltzmann collision 
velocity distribution. Its variation with temperature is 
displayed in Fig. [3] for the RA of with para-H2 
{vq — 0,jo — 0) and with ortho-Il2 {vq — 0, jo = !)■ As a 
follow-up of Fig. [21 the rate is found about 4 times larger 
in the ortho-H2 case than in the para-H2 case. Note that 
the resonances in the latter case are too narrow to in- 
fluence the rate above 200 K. The relative abundances 
of ortho-H2 and para-Il2 could be out of thermal equi- 
librium in the interstellar clouds The obtained RA 
rates suggest that the ratio of ortho-H2— H~ to para- 
H2— H~" is enhanced by a factor of 4 compared to the 
ratio of ortho-H2 to para-H2. 
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FIG. 3. The rate coefficient Ara as a function of temperature 
for RA of H~ with para-H2 (wo ~ 0,jo = 0) (fuU hne) and 
with ortho-H2 {vo = 0,jo = 1) (broken hne). 



V. DISCUSSION: IS THE FORMATION OF H3 
IN THE ISM POSSIBLE? 

The production rate of H3 in the ISM is directly re- 
lated to the presence of H~. Here, we propose a rough 
estimation of the abundance of H" in the ISM assum- 
ing that it is formed only by dissociative attachment of 
an electron to H2. Other possible mechanisms such as 
radiative attachment to H are neglected. The chemistry 
of interstellar clouds is initiated by ionization of molecu- 
lar hydrogen by cosmic rays with a typical rate constant 

~ 3 X 10~^^s~^ in diffuse interstellar clouds [1]. (Cloud 
densities are ~ 10^ cm""^ in diffuse clouds and ^ 10'' 



cm~ in dense clouds.) The ionized molecular hydrogen 
quickly forms in collisions with H2, with a rate 
constant ^ 2 x 10~^cm'^/s The electron escaped after 
ionization of H2 has a large kinetic energy and undergoes 
many elastic collisions with environmental II2 before it 
thermalizes. In each elastic collision with II2, the electron 
looses a fraction (about Ame/mn^) of its incident energy. 
For example, the electron should experience about 6200 
elastic collisions with H2 before its energy is decreased 
from 10 eV to 10 meV (T ~120 K). Possible inelas- 
tic e^+Il2 collisions will lead to vibrational excitation 
of II2 and to dissociative attachment (DA), e~+ll2 — 
H-|-H~. The DA reaction is allowed for collision ener- 
gies above the threshold at 3.7 eV with the cross section 
about aoA 3 x lO'^i cm^ at 4 eV [Ullll. 

We can roughly estimate the fraction /da of escaped 
electrons that would form H~ rather than undergo the 
thermalization process. This gives us an estimate about 
the production rate of H~ in the ISM. The upper bound 
for the e~+ii2 elastic cross-section a^i is about 1.7 x 
10-15 cm^ [sil. If we take Nth = 1000 elastic colli- 
sions corresponding to the thermalization down to the 
energy below which DA is impossible, we should com- 
pare a el I Nth with a DA- It gives approximately the value 
/da = c^DANth/ c^ei ^ 0.0018 for the fraction of escaped 
electrons that form H" by dissociative attachment to H2 . 
Therefore, the rate of H" production in cm"^ in the ISM 



can be estimated as C"(H2)/da = Cn(H2)<^DANth/<Tei- 
The H" ion has only one bound electronic state and can- 
not be detected directly. However, if it forms a molecular 
ion AH~ by an association with an atom or molecule A, 
its presence in the ISM can be proven indirectly by rovi- 
brational absorption spectroscopy of AH~ . 

One of the motivation for this study was the investi- 
gation whether H;3' can be formed in the ISM and be 
detected by photoabsorption spectroscopy. The stability 
of Hj^ was confirmed in this study. Therefore, Hj^ can 
exist in the ISM in cold clouds (T < 100 K). However, to 
be observed by far-infrared absorption spectroscopy, H;^ 
should be present in the ISM in relatively large amounts. 
The Hg molecular ion and the H~ ion can be destroyed 
in the ISM by cosmic rays and by mutual neutralization 
with positive ions. These processes limit the absolute 



abundance of H 



3 



In order to determine the H3 abun- 



dance, one has to consider rate equations for all reac- 
tions involving formation and removal of H~ and H3" in 
the ISM. However, a rough estimation about H^ and H3 
abundance in the ISM can be made using the available 
data. 

First, we estimate the H~ abundance n(H~). Accord- 
ing to our model the rate of H~ production in 1 cm'^ is 
given by C"(H2)/da- The principal channel of destruc- 
tion is probably due to the mutual neutralization with 
positive ions. The number density of positive ions can 
be taken to be equal to the number density of electrons in 
the ISM. In diffuse clouds, the electron number density 
is about 0.1 % of n(H2), i.e. n+ ~ 0.01 cm"'^. Therefore, 
the rate of removal of H~ from the ISM is fc±n(H~)n_(- (in 
1 cm'^), where k± is the rate coefficient for the mutual 
neutralization. We take the value of k+ ^ 10~^cm'^/s 
for the H--hH+ ^H-^H reaction jH-fsil at 10 meV. 
Therefore, we derive the equilibrium abundance of H^: 



n(H-) = C"(H2)/DA/(fc±n+) ^ 5 x 10" 



If we 



take the size of a cloud to be 10 pc, the resulting column 
density is ~ 10^'^ cm~^, which is a reasonable number for 
observation. 

The estimation of the H3" abundance n{}i'^) can be 
made in a similar way. The rate of H3" removal is 
determined by the similar formula: fc±n(H^)rt+. The 
rate of H3" formation is A:RA'T'(H~)n(H2). The equi- 
librium H^ abundance is then given by ^(Hg") = 
/cRA'T-(II~)n(H2)/(fc±n+) ^ 10~i^cm~'^. Such an abun- 
dance would produce a column density, which is small 
even for a large interstellar cloud. 

If we combine the formulas for H~ and H3" abundances 
we obtain 



fcRAC/RA f n{ll 
(fc±)2 



(37) 



The formula suggests that the Hg" abundance should be 
larger in an environment with smaller degree of ionization 
(i.e. with a smaller number density of positive ions). 
The environment should be cold enough (<100 K) for 
the H;^ ion to be stable with respect to collisions with 
other species. Therefore, if H;^ can be detected in the 
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ISM, one has to search for it in cold dense intersteUar 
clouds. If is detected it would be a proof that is 
also present in the ISM. 
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